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Gauge-invariant spatial correlations between two and three quarks inside a hadron are measured within 
quenched and unquenched QCD. These correlators provide information on the shape and multipole moments 
of the pion, the rho, the nucleon and the A. 

1. Introduction 



Gauge-invariant two- and three-density corre- 
lators inside a hadron reduce to the square of the 
wave function in the non-relativistic limit, yield- 
ing detailed information on hadron structure. 
Quark spatial distributions, hadronic shapes, 
charge radii, etc., can be extracted. An interest- 
ing question is whether the nucleon is deformed. 
Strong evidence for deformation in the nucleon 
and/or A is provided by recent accurate mea- 
surements ^ in photoprodution experiments on 
the nucleon. One can compare with the experi- 
mental results by calculating the N to A"*" tran- 
sition form factors |2]. Moreover, the A+ defor- 
mation can be directly obtained from the wave 
function. The nucleon, however, has zero spec- 
troscopic quadrupole moment, since it is a spin 
1/2 particle. This implies that any deformation 
that intrinsically may be present for a given back- 
ground gauge field will average out to zero in the 
gauge configuration ensemble. 

2. Observables 

In this work we consider equal-time correlators. 
For a meson we calculate the matrix element 

Cr(r,i)- J dV (M|p«(r',t)p'^(r' + r,t)|M)(l) 

with Pr(r,i) u{r,t)Tu{r,t) :. We take F = 
7o for the normal ordered charge density oper- 
ator, and 1 for the matter density. In the non- 
relativistic limit this allows to extract respectively 



the charge and matter density distributions. Un- 
hke Bethe-Salpeter amplitudes, these correlators 
are gauge- invariant. Here we take the u— and d— 
quarks to be degenerate in mass. 

The form factors F{q^) at low momenta can 
be extracted from the Fourier transform of the 
charge density-density correlator 



d'^r exp (iq.r) C^g (r, t) /2m^ 



*Talk presented by C. Alexandrou. 



^ . (2) 

For baryons the charge distribution is obtained 
by using three density insertions. This involves 
two relative distances and it is computed effi- 
ciently by using EFT. The diagrams involved are 
shown in Fig. |21 The top right diagram requires 
the all-to-all propagator 4 and will not be con- 
sidered here. We check that the contribution of 
this diagram is small by comparing the integrated 
one-particle density / d^r2C{ri,r2,t) with the 
two-density correlator, / d'^r' {h\p'^{r' ,t)p^{r' + 
ri,t)\h), shown by the lower diagram of Fig. |21 

3. Results 

We first discuss the results obtained in the 
quenched approximation at /3 = 6.0 for a lattice of 
size 16'^ X 32. The parameters of the simulation 
for both quenched and unquenched calculations 
are summarized in the following Table. 

The charge and matter density distributions for 
the pion and the rho mesons arc shown in Fig.[21at 
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Figure 1. Three-density correlator for a baryon 
where all the insertions are taken at equal times 
t. t and T/2 — t are taken large enough to filter 
out the excited baryonic states. 

Quenched Unquenched(from pi) 

K rriTr/mp #confs. Ksea m-nlrrip #confs. 

0.153 081 220 (U56 083 150"" 

0.154 0.78 220 0.157 0.76 200 



K = 0.154. We find that the charge distribution 
is broader than the matter distribution, in agree- 
ment with Ref. 'SI. The matter distributions for 
the pion and the rho are almost identical, but the 
density distribution is considerably broader for 
the rho. As seen in Fig. |31 the nucleon and the 
A"*" have very similar matter distributions to the 
rho meson, whereas larger variations are visible 
in the charge density distributions. 

In Fig. 21 we show the pion form factor which, 
being the simplest, can be reliably extracted from 
our set of data. Its low momentum behaviour 
is well described by the vector dominance re- 
sult F(q^) ^ 1/(1-1- c^/m?). We observe a 
weak dependence on k obtaining 0.57, 0.59 and 
0.61 GeV at K = 0.153,0.1554 and 0.155 respec- 
tively. Extrapolating linearly in we obtain 
m ^ 0.65 GeV within 15% of the rho meson mass. 

Hadron deformation can be investigated by 
evaluating the root mean squared (rms) radius 
along the spin axis and perpendicular to it. We 
show the charge rms radii in Fig. as a function 
of the quark mass. The pion, having equal rms 
radii, is spherical whereas for the rho the asymme- 
try between the longitudinal and transverse radii 
increases as we approach the chiral limit. A reli- 
able extrapolation to the chiral limit needs lighter 
quarks and a larger lattice to contain the rho. A 
similar analysis for the nucleon shows that it is 
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Figure 2. Charge and matter density distribu- 
tions at K = 0.154 for the pion and the rho. 
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Figure 3. Charge and matter density distribu- 
tions at K = 0.153. 

spherical as expected, whereas the A+ shows no 
statistically significant deformation W. The rms 
radii computed from the matter distribution show 
no such deformation also in the case of the rho. 

We now compare our quenched results for 
the charge correlator with those obtained us- 
ing two degenerate flavors of Wilson fermions, 
with similar lattice spacing and rriT^/mp values. 
Fig. shows that unquenching leaves the pion 
unchanged, whereas the rho, the nucleon and the 
A"*" become broader. The asymmetry in the rho 
increases as seen by comparing the longitudinal 
and transverse radii in Fig. [3 The A+ shows a 
small asymmetry which points to pion cloud con- 
tributions to the deformation. 

The three-density charge correlator was com- 
puted in the quenched approximation using 30 
configurations at k = 0.15 and 0.154. More de- 
tailed information on hadron structure can be ex- 
tracted from this correlator. As an example we 
show in Fig. |7| the u- and d- quark spatial dis- 
tributions. The systematically broader d-quark 
distribution gives rise to a negative neutron rms 
charge radius. 
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Figure 4. 1/F(q^) versus at k = 0.154. The 
solid line is a fit to the form (1 + /m?). 
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Figure 5. (z^) and {{x^ +y'^)/2) vs the pion mass 
squared. The lines are linear fits to the quenched 
data. The empty circles show full QCD results. 

4. Conclusions 

We have presented a gauge-invariant determi- 
nation of hadron charge and matter density dis- 
tributions. The main conclusions from this study 
are the following: 1) The matter density distri- 
bution is very similar for the tt, p, nucleon and 
A"*", unlike the charge density which is narrower 
for the pion than for the rest. In all cases the 
charge density distribution is broader than the 
matter density distribution. 2) Vector dominance 
provides a good description of the pion form fac- 
tor at low momenta. 3) Unquenching leaves the 
pion size unchanged but the rho, the nucleon and 
A+ become broader. 4) In the quenched approxi- 
mation the rho in the 0-spin projection is prolate 
whereas the A+ has no statistically significant de- 
formation. 5) Unquenching increases the rho de- 
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Figure 6. (a) Density-density correlators, C-yo(r), 
for the pion, the rho, the nucleon and the A"*" at 
K = 0.154 vs |r|. (b) Same as (a) but with two 
dynamical quark flavors at k = 0.157. Errors bars 
are omitted for clarity. 
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Figure 7. The u— and d— distributions inside the 
neutron for k ~ 0.15 and 0.154. 

formation and produces a small deformation for 
the A+ with the +3/2 spin state being slightly 
oblate. 6) The d-quark spatial distribution in the 
neutron extracted from the three density corre- 
lator is broader than that of the u-quark, thus 
accounting for the negative rms charge radius of 
the neutron. 
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